QCD HARD SCATTERING RESULTS FROM PHENIX 
AT RHIC 



David d'Enterria for the PHENIX Collaboration 

Nevis Laboratories, Columbia University 
Irvington, NY 10533, and New York, NY 10027, USA 
denterria@nevis.columbia.edu 



Abstract Data on hadron production at high transverse momentum (p-r > 2 GeV/c) in 
p+p, d+Au, and Au+Au collisions at v /s^iv = 200 GeV from the PHENIX ex- 
periment at RHIC are reviewed. The single inclusive spectrum of light hadrons 
produced in central Au+Au reactions shows significant differences compared to 
p+p and d+Au collisions, and provides interesting information on the proper- 
ties of the underlying QCD medium present in heavy-ion reactions at collider 
energies. 
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1. Introduction 

The fundamental degrees of freedom of the theory of the strong interaction 
(Quantum Chromodynamics, QCD) are colored quarks (fermions) and gluons 
(gauge bosons). In the real world, however, quarks and gluons are not observed 
free but confined in color singlet states (hadrons). In addition, the observed ef- 
fective constituent mass of the quarks (fh Ui d ~ 300 MeV for the lightest u and 
d quarks) is much larger than the bare current mass of the QCD Lagrangian 
( m u,d ~ 2 — 8 MeV). Theoretically understanding from first principles these 
two basic properties of QCD: confinement and chiral symmetry breaking, re- 
mains in fact one of the most important problems in fundamental physics [1]. 
Lattice calculations of QCD in the bulk [2] predict that above an energy density 
of e cr it ~ 0.7 ± 0.3 GeV/fm 3 , the (long range part of the) strong color poten- 
tial becomes screened and, as a result, quarks and gluons can move freely with 
their bare mass within the deconfined medium. Experimentally accessing and 
studying the properties of this "Quark-Gluon Plasma" (QGP) phase is one of 
the main goals of the physics program of the PHENIX experiment at the Rela- 
tivistic Heavy-Ion Collider (RHIC) at BNL. By colliding heavy nuclei at high 
energies one expects to produce a thermalized system of quarks and gluons 
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with energy densities above e cr it, albeit for very short time scales, O(10 23 
s) and small volumes, O(10 2 fm 3 ). In the period 2000-2003, PHENIX col- 
lected data from Au+Au, d+Au and p+p collisions at y/s NN = (130)200 GeV. 
The most interesting phenomena observed so far in the data are in the high px 
(Pt 2 GeV/c) sector where the production of hadrons in (non-peripheral) 
Au+Au collisions shows substantial differences compared to more elementary 
reactions either in free space (p+p, e + e~) or in a "cold" nuclear matter en- 
vironment (d+Au). Hard scattering processes are an excellent experimental 
probe in heavy-ion collisions inasmuch as: (i) they are the result of violent 
(large Q 2 ) short-distance (parton-parton) collisions in the first instants of the 
reaction (r ~ 1/pt < 0.1 fm/c) and, as such, direct probes of the partonic 
phase(s) of the reaction; (ii) in the absence of medium effects, their cross sec- 
tion in A+A reactions is expected to scale simply with that measured in p+p 
collisions times the number of scattering centers (A 2 , see discussion below); 
and (iii) their production yields can be theoretically calculable in QCD via 
standard perturbative (collinear factorization plus Glauber multiscattering [3]) 
or via "classical-field" ("Color-Glass-Condensate" approach [4]) methods. 

2. The PHENIX experiment 

The PHENIX detector [5] is specifically designed to measure hard QCD 
probes such as high px hadrons, direct photon radiation, lepton pairs, and 
heavy flavour production. PHENIX achieves good mass and PID resolution, 
and small granularity by combining 13 detector subsystems (~350,000 chan- 
nels) divided into: (i) 2 central arm spectrometers for electron, photon and 
hadron measurement at mid-rapidity (\r}\ < 0.38, A</> = vr/2); (ii) 2 forward- 
backward (\r]\ = 1.15 - 2.25, A(f) = 2tt) spectrometers for muon detection; 
and (iii) 4 global (inner) detectors for trigger and centrality selection. Neu- 
tral mesons are reconstructed through invariant mass analysis of decay 7 pairs 
detected in two types of electromagnetic calorimeters (15552 lead scintillator 
towers with ISXq, and 9216 lead glass modules with 14.4Xo). The trajectories 
and momenta of charged hadrons in the axial central magnetic field (B max = 
1.15 T m) are measured by a drift chamber (DC) and three layers of MWPC's 
with pad readout (PC). Hadron identification (tt^, K^, and p,p) is achieved 
by matching the reconstructed tracks to hits in a time-of-flight wall (TOF). 

The occurrence of a p+p, d+Au or Au+Au collision (with vertex position 
\z\ < 60 cm within the center of the detector) is triggered by a coincidence be- 
tween the two Beam-Beam Counters (BBC) 1 which cover \rj\ = 3.0 - 3.9. These 
minimum bias triggers accept respectively (52±4)%, (88±4)%, and (92±3)% 
of the total inelastic cross-sections (cr™ el ~ 42 mb, (f^ u « 2150 mb, and 

'Plus the Zero Degree Calorimeters, ZDC, at |0| < 2 mrad in the case of Au+Au. 
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° % AuAu ~ 6850 mb). Since hard probes are rare, high luminosities are required. 
PHENIX is an efficient high interaction rate experiment with a state-of-the-art 
data acquisition system capable of recording 120 MB/s to disk with event sizes 
of ~100 KB, and event rates of ~l-2 KHz. A total of 0.2-10 9 , 5.5-10 9 , and 
4.0- 10 9 events have been sampled in Au+Au, d+Au, and p+p collisions respec- 
tively at ^/s^ = 200 GeV. 

3. Scaling of yields and cross-sections from p+p to A+B 

Insights on the mechanisms of particle production (and "destruction") in 
nucleus-nucleus (A+B) collisions are obtained from the study of the scaling 
behavior of their yields with respect to p+p collisions. Depending on the px 
range in which they are produced, two different "scaling laws" are relevant: 

• Soft particle production at low pr (pr ^5 1 GeV/c) in A+B reactions is 
dominated by non-perturbative nucleon-nucleon collisions with small momen- 
tum transfers. Since those collisions occur in "parallel" and with large proba- 
bilities, they suffer destructive interference which effectively limits the maxi- 
mal number of processes that can lead to secondary particle production [6]. As 
a result soft particle yields turn out to be proportional to the average number 
of participating (or "wounded") nucleons, i.e. nucleons which undergo at least 
one inelastic collision [7]. Thus, in a A+B reaction at impact parameter b: 

where N part {b) can be obtained from the nuclear thickness functions T^ ; b(6) = 
J dz Pa,b(z, b) (normalized to A, B resp.) via a Glauber model 2 . Such a "par- 
ticipant scaling" is indeed approximately observed in the total multiplicities 
measured in heavy-ion collisions [8]. Integrating N part (b) over impact param- 
eter one obtains the average number of participant nucleons in a minimum bias 
A+B collision: (N part }\ MB = (Aa l p n £ l +B <r p n A el ) / cr 1 ^ 1 . Using this expression, 

Eq. (1), and dN AA = da AA /a l AA , the soft invariant yields and cross-sections 
in a symmetric A+A minimum bias collision 3 are related to the corresponding 
p+p ones via: 
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• Hard particle production at high p? (pr ^ 2 GeV/c) results from incoher- 
ent parton-parton scatterings with large Q 2 . In this regime, the pQCD "fac- 



2 The corresponding formulas can be found e.g. in [7]. 

3 Also, for min. bias A+A collisions: (ATj, art )| MB = 2Aa 1 p n J' 1 ja^f fs 2AttR a /tt(2R a ) 2 A/2. 
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torization theorem" [9] holds and the inelastic cross-section for the produc- 
tion of a given particle can be separated in the product of long-distance (non- 
perturbative partem distribution functions, f a /A, and fragmentation functions, 
D c / h ) and short-distance (parton-parton scattering) contributions: 



j„hard 
jji_AB—>h 



fa/A(x,Q 2 )® fb/B(.X,Q 2 ) 



j hard 
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D c/h (z,Q 2 ). (3) 



d3p ,~, ,*v~,<* ,~ d3p 

The assumption of incoherent scattering at high pr entails also that f a /A = 
A - f a / N , i.e. that the density of partons in a nucleus A should be equivalent to 
the superposition of A independent nucleons. Thus, 

da h A a X± 
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and minimum bias hard cross-sections in A+B are expected to scale simply as 

(5) 
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In the general case, for a given A+B reaction with impact parameter b: 
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d 3 p 

where TAB{b) = f d 2 sTA(s)TB(b — s) (normalized to AB) is the nuclear 
overlap function 4 at b. Since Tab is proportional to the number of nucleon- 
nucleon (NN) collisions: Tab (6) = N co u(b) / a™ el , one alternatively quotes 
Eq. (6) in the form of "(binary) collision scaling" of invariant yields: 
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From (5) and (7), it is easy to see that 5 (iV co ^) | MB = AB 

Following Eqs. (2) and (5), since (J % v n A d AA oc R 2 A oc A 2 / 3 , the atomic number 
dependence of soft processes in symmetric A+A reactions is of the type 
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whereas the A dependence of hard processes is of the form: 
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4 The "natural" magnitudes of T A (b) and T AB (b) are A/irR A and AB/n(R A + R B ) 2 respectively [10]. 
5 Also, for min. bias A+A collisions: {N co11 )\ M b = A 2 o™ 1 1 o^f Ri A 2 a™ el /n(2R A ) 2 A 4 / 3 /4. 
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4. High p T production: Au+Au vs. p+p 

One of the most interesting experimental results at RHIC so far is the break- 
down of the expected incoherent parton scattering assumption for high pt pro- 
duction in non-peripheral Au+Au collisions. Fig. 1 shows the comparison 
of the measured p+p 7r° spectrum [11] to peripheral (left) and central (right) 
Au+Au spectra [12], and to standard NLO pQCD calculations [13]. Whereas 
peripheral data is consistent with a simple superposition of individual NN 
collisions, central data shows a suppression factor of 4 - 5 with respect to this 
expectation. 
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Figure 1. Invariant tv° yields measured by PHENIX in peripheral (left) and in central (right) 
Au+Au collisions (stars) [12], compared to the Tab scaled p+p 7r° cross-section (circles) [11] 
and to a NLO pQCD calculation (gray line) [13]. The yellow band around the scaled p+p points 
represents the overall normalization uncertainties. 

From Eq. (6), it is customary to quantify the medium effects at high pt using 
the nuclear modification factor given by the ratio of the A+A invariant yields 
to the p+p cross-sections scaled by the nuclear overlap at impact parameter b: 

d^Nf A (b)/dydp T 

Kaa{PT,0) - j= 7T77 2 71 , • (10) 

(T AB (b)) ■ d 2 a^ p /dydpT 

Raa{pt) measures the deviation of A+A from an incoherent superposition of 
NN collisions in terms of suppression (Raa <1) or enhancement (Raa >!)• 
Figure 2 shows Raa as a function of pt for several tt° measurements in high- 
energy A+A collisions. Much of the excitement at RHIC comes from the fact 
that the PHENIX Raa values for central collisions at 200 GeV (circles) and 
130 GeV (triangles) are noticeably below unity in contrast to the enhanced pro- 
duction (Raa >1) observed at CERN-ISR [14] (stars) and CERN-SPS [15] 
(squares) energies. This enhanced production, observed first in p+A fixed- 
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target experiments [16] ("Cronin effect"), is interpreted in terms of multiple 
initial-state soft and semi-hard interactions which broaden the transverse mo- 
mentum of the colliding partons prior to the hard scattering itself. 
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Figure 2. Nuclear modification factor, Raa (pt), for 7r° measured in central ion-ion reactions 
at CERN-SPS [15], CERN-ISR [14], and BNL-RHIC [17, 12] energies. The dashed (dotted) 
line is the expectation of "N co ii (N par t) scaling" for "hard" ("soft") particle production. 

The breakdown of the expectations from collinear factorization for high pj> 
production in central A+A collisions at RHIC by such a large factor, has been 
interpreted as due to one (or more) of the following facts: 

1 Breakdown of leading-twist QCD collinear factorization itself: The in- 
coherence between long- and short-distance effects in which the factor- 
ized product Eq. (3) relies upon, does not hold for A+A collisions. 

2 Strong initial-state effects: The parton distribution functions in the nu- 
clei are strongly modified: f a i A « A ■ f a / p in the relevant (x,Q 2 ) 
range, resulting in a reduced number of effective partonic scattering cen- 
ters in the initial-state. 

3 Strong final-state effects: The parton fragmentation functions (or, more 
generally, any post hard collision effect on the scattered partons) are 
strongly modified in the nuclear medium compared to free space. 

Explanations 1. and 2. are usually invoked in the context of the "Color- 
Glass-Condensate" picture [4] which assumes that the kinematical conditions 
prevailing in the initial-state of an atomic nucleus boosted to RHIC energies 
(moderate Q 2 rs Q 2 = 1-2 GeV 2 /c 2 , and small parton fractional momenta 
x < 1CP 2 ) are such that nonlinear QCD effects (g + g — > g processes, am- 
plified by a A 1 / 3 factor compared to the proton case) are important and lead 
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to a strong saturation of the partem (mostly gluon) densities in the nuclei. In 
this scenario, one expects N part (instead of N co u) scaling at moderately high 
Pt's [18], as approximately observed in the data (dotted line in Fig. 2). Ex- 
planation 3., on the other hand, relies on the expectations of "jet quenching" 
in a Quark Gluon Plasma [19] in which the hard scattered partons lose en- 
ergy by final-state "gluonstrahlung" in the dense partonic system formed in the 
reaction. After traversing the medium, the partons fragment into high p-r (lead- 
ing) hadrons with a reduced energy compared to standard fragmentation in the 
"vacuum". Different jet quenching calculations can reproduce the magnitude 
of the 7T° suppression assuming the formation of a hot and dense partonic sys- 
tem characterized by different, but closely related, properties: i) large initial 
gluon densities dN 9 jdy ps 1000 [20], ii) large "transport coefficients" % ps 3.5 
GeV/fm 2 [21], iii) high opacities L/X ps 3.5 [22], iv) effective partem energy 
losses of the order of dE/dx ps 14 GeV/fm [23], or v) plasma temperatures of 
T ps 0.4 GeV [24]. 

5. High p T production: d+Au vs p+p 

In order to disentangle between the two different (QGP and CGC) QCD 
scenarios, it was of paramount interest to determine experimentally the modi- 
fication, if any, of high px hadron production due to initial-state nuclear effects 
alone, i.e. for a system in which a hot and dense medium is not produced in 
the final state. The results of high pt vr° production at midrapidity in d+Au 
collisions at ^/s NN = 200 GeV [25] do not show any indication of suppres- 
sion (Fig. 3). On the contrary, 7r° production seems to be slightly enhanced 
(RdAu ~ 1-1) compared to the expectations of collinear factorization. 
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Figure 3. Nuclear modification factor, Raa(pt), for it measured at y = in minimum bias 
d+Au (squares) [25] and central Au+Au (circles) [12] reactions at ^/s NN = 200 GeV. 
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This result, reminiscent of the "Cronin enhancement", indicates that the 
observed suppression in Au+Au central collisions at mid-rapidity is not an 
initial-state effect arising from strong modifications of the gluon distribution 
functions in nuclei at moderately small values of partem fractional momenta x, 
but more likely a final-state effect of the produced dense medium. 

6. Summary 

PHENIX has measured 7r° at mid-rapidity up to pr ps 10 GeV/c in p+p, 
d+Au and Au+Au collisions at ^/s~^ = 200 GeV. The spectral shape and in- 
variant yields in peripheral Au+Au reactions are consistent with those of p+p 
reactions scaled by the number of inelastic NN collisions, in agreement with 
pQCD collinear factorization expectations. In central Au+Au reactions, the 
7T° are suppressed by a factor 4-5 with respect to the same expectations. The 
magnitude of this deficit can be reproduced by pQCD-based partem final-state 
energy loss calculations in an opaque medium, as well as by initial-state gluon 
saturation models. The unquenched high pr vr° production in minimum bias 
d+Au collisions, however, seems to rule out any significant "cold" nuclear 
matter effect as responsible for the observed suppression in Au+Au central 
reactions. 
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